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Charm-Quark Contribution to Kl — > fJ, + n~ at Next-to-Next-to-Leading Order 
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We calculate the charm-quark contribution to the decay Kl —* ^ + p-~ m next-to-next-to-leading 
order of QCD. This new contribution reduces the theoretical uncertainty in the relevant parameter P c 
from ±22% down to ±7%, corresponding to scale uncertainties of ±3% and ±6% in the short-distance 
part of the branching ratio and the determination of the Wolfenstein parameter p from Kl — > p, + p~ ■ 
The error in P c = 0.115±0.018 is now in equal shares due to the combined scale uncertainties and the 
current uncertainty in the charm-quark mass. We find B(Kl ~* /i + /i~)sD = (0.79 ±0.12) x 10~ 9 , 
with the present uncertainty in the Cabibbo-Kobayashi-Maskawa matrix element Vtd being the 
dominant individual source in the quoted error. 

PACS numbers: 12.15.Hh, 12.38.Bx, 13.20.Eb 



The study of the rare process Kl — > p + p~ has played 
a central role in unraveling the flavor content and struc- 
ture of the standard model (SM) of electroweak interac- 
tions 0. These glory days have passed, but still today 
Kl — * P + P~ provides useful information on the short- 
distance dynamics of |AS| = 1 flavor-changing-neutral- 
current transitions despite the fact that its decay ampli- 
tude is dominated by the long-distance two photon con- 
tribution Kl — > 7*7* — > P + P- While the absorptive 
part of the latter correction is calculable with high preci- 
sion in terms of the Kl — ► 77 rate the corresponding dis- 
persive part represents a significant source of theoretical 
uncertainty. In fact long- and short-distance dispersive 
pieces cancel against each other in large parts and the 



J 



measured total Kl — > A* + A* _ ra ^ e \M I s nearly saturated 
by the absorptive two photon contribution. The precision 
in the determination of the dispersive pieces therefore 
controls the accuracy of possible bounds on the real part 
of the Cabibbo-Kobayashi-Maskawa (CKM) element Vtd 
or, equivalcntly, the Wolfenstein parameter p. In view 
of the recent experimental Q and theoretical Q devel- 
opments concerning the dispersive long-distance part of 
the Kl — > p + p~ decay amplitude it is also worthwhile to 
improve the theoretical accuracy of the associated short- 
distance contribution. This is the purpose of this Letter. 

The branching ratio for the dispersive short-distance 
part of Kl — > p + p~ can b e written as Q 



B(K L -> /i+aOsd 
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7T Z sm 



t(K+) 
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0.225 
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where Aj = V* s Vid denote the relevant CKM factors. 
There is also a short-distance two-loop electroweak con- 
tribution in the two-photon mediated decay amplitude 
0. Following 0, where this contribution is included 
in the two-photon correction itself, we do not add it to 
the short-distance contribution in Eq. ft}. The apparent 
strong dependence of B(Kl — * A* + aOsd on A = |K»s| 
is spurious as P c is proportional to 1/A . In quoting 
the value for P c we will set A = 0.225. The electro- 
magnetic coupling a and the weak mixing angle sin 2 9 W 
entering B(K L — > M + ' p~) are naturally evaluated at the 
electroweak scale [7j| . Then the leading term in the heavy 
top expansion of the electroweak two-loop corrections to 
Y(xt) amounts to typically —1.5% for the modified mini- 



mal subtraction scheme (MS) definition of a and Sm U\\r 
[8j • In obtaining the numerical value of Eq. J2J) we have 
employed a = a m (M z ) = 1/127.9, sin 2 6 W = sin 2 9^ = 
0.231, and B{K+ -> p+v^) = (63.39 ± 0.18) x 10" 2 @. 

The function Y(xt) in Eq. depends on the top quark 
MS mass through xt = mf(pt)/M^. It originates from 
Z-penguin and electroweak box diagrams with an inter- 
nal top quark. As the relevant operator has a vanish- 
ing anomalous dimension and the energy scales involved 
are of the order of the electroweak scale or higher, the 
function Y(xt) can be calculated within ordinary per- 
turbation theory. It is known through next-to-leading 
order (NLO) [13, EH, 

with a scale uncertainty due to 
the top quark matching scale pt — 0(m t ) of slightly 
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TABLE I: Input parameters used in the numerical analysis of 

P c , B(K L -> ^ + ai _ )sd, and p. 
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Parameter 



Value ± Error 



Reference 



m c (m c ) [GeV] 
a s {M z ) 



1.30 ±0.05 
0.1187 ±0.0020 



ReA t 
ReAc 



[io- 



-3.11 



+0.13 
-0.14 



-0.22098i°;S 



QM, our average 

m 

[15] 



more than ±2%. Converting; the top quark pole mass 
of M t = (172.5 ± 2.3) GeV [lj at three loops to m t {M t ) 
[13 and relating m t (M t ) to m t {m t ) = (162.8 ± 2.2) GeV 
using the one-loop renormalization group (RG), we find 
Y(xt) = 0.950 ± 0.049. The given uncertainty combines 
linearly an error of ±0.029 due to the error of mt(mt) and 
an error of ±0.020 obtained by varying /i t in the range 
60 GeV < n t < 240 GeV. 

The calculable parameter P c entering Eq. results 
from Z-penguin and electroweak box diagrams involv- 
ing internal charm-quark exchange. As now both high- 
and low-energy scales, namely, [i w = 0(M W ) and /i c = 
0(m c ), are involved, a complete RG analysis of this term 
is required. In this manner, large logarithms hiQi^/fi^) 
are resummed to all orders in a s . The large scale uncer- 
tainty due to n c of ±44% in the leading order result was 
a strongmotivation for the NLO analysis of this contri- 
bution BEl 

Performing the RG running from fj, w down to /i/, = 
0(rrib) in an effective five-flavor theory and the subse- 
quent evolution from ^ down to /j, c in an effective four- 
flavor theory, we obtain at NLO 



... i .... i .... i .... i .... i .... i .... i .. . 

1 1.25 1.5 1.75 2 2.25 2.5 2.75 3 

A*c [GeV] 



(0.106 ±0.034) 



0.225 
A 



(3) 



where the parametric errors correspond to the ranges of 
the charm-quark MS mass m c (m c ) and the strong cou- 
pling constant a s (M z ) given in Tablcfl] The final er- 
ror has been obtained by performing a detailed analy- 
sis of the individual sources of uncertainty entering the 
NLO prediction of P c using a modified version of the 
CKMFITTER package [15| . The same statistical treatment 
of errors will be applied in Eqs. (J2J, and 

The dependence of P c on \i c can be seen in Fig.^ 
The solid line in the upper plot shows the NLO result 
obtained by evaluating a s (p c ) from a s (M z ) solving the 
RG equation of a s numerically, while the dashed and 
dotted lines are obtained by first determining the scale 
parameter A^jg from a,(M z ), either using the explicit 
solution of the RG equation of a s or by solving the RG 
equation of a s iteratively for A^g, and subsequently cal- 
culating a s {n c ) from A^jg. The corresponding two- loop 
values for a s (p c ) have been obtained with the program 




1.75 2 2.25 

A*c [GeV] 

FIG. 1: P c as a function of p c at NLO (upper plot) and NNLO 
(lower plot). The three different lines correspond to three differ- 
ent methods of computing a s (/i c ) from a 3 (M z ) (see text). 



RUNDEC ^(| . Obviously, the difference between the three 
curves is due to higher order terms and has to be re- 
garded as part of the theoretical error. With its size of 
±0.006 it is almost comparable to the variation of the 
NLO result due to fi c , amounting to ±0.016. In |5j a 
larger value for the latter uncertainty has been quoted. 
The observed difference is related to the definition of the 
charm-quark mass. Replacing m c (m c ) in the logarithms 
ln(/x^/m^) of the one-loop matrix elements by the more 
appropriate m c (fx c ) leads to a significant reduction of the 
dependence of P c on u r . A detailed discussion of this is- 
sue can be found in [l7j. Finally, while in only /jl c 
was varied, the theoretical error given in Eq. © includes 
also the dependence on /it> and fi w of combined ±0.001. 
The specified scale uncertainties correspond to the ranges 
1 GeV < \l c < 3 GeV, 2.5 GeV < [i h < 10 GeV, and 
40 GeV < n„ < 160 GeV. 

Using the input parameters listed in Tabled we find 
from Eqs. ©-© at NLO 



B{K L -> a* + aOsd = (0.77 ± 0.08 Pc ± 0.08, 
= (0.77 ±0.16) x 10~ 9 



'other) X 10 9 

(4) 



where the second error in the first line collects the uncer- 
tainties due to k m , Y(xt), and the CKM elements. 
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FIG. 2: Examples of Feynman diagrams arising in the full 
SM (left column), describing the mixing of operators (cen- 
ter column) and the matrix elements (right column) in the 
Z-penguin (upper row) and the electroweak box (lower row) 
sector. Only the divergent pieces of the diagrams displayed 
in the center column have to be computed, while the Feyn- 
man graphs shown on the left- and right-hand side are needed 
including their finite parts. 



As the uncertainties in Eqs. J2J) and 0} coming from 
M t , m c (m c ) and the CKM parameters should be de- 
creased in the coming years it is also desirable to reduce 
the theoretical uncertainty in P c . To this end, we here 
extend the NLO analysis of P c presented in jfl to 
the next-to-next-to-leading order (NNLO). This requires 
the computation of three-loop anomalous dimensions of 
certain operators and of certain two-loop contributions. 

The main components of the NNLO calculation, which 
aims at resumming all 0(a™ \n n ~ 1 (p 2 v /fj,D) logarithms in 
P c , are (i) the 0(a 2 s ) matching corrections to the relevant 
Wilson coefficients arising at fi W) (ii) the 0(a\) anoma- 
lous dimensions describing the mixing of the dimension- 
six and -eight operators, (Hi) the 0(a 2 ) threshold cor- 
rections to the Wilson coefficients originating at /if,, and 
(iv) the 0(a 2 s ) matrix elements of some of the operators 
emerging at /i c . To determine the contributions of type 
(i), (Hi), and (iv) one must calculate two-loop Green 
functions in the full SM and in effective theories with 
five or four flavors. Sample diagrams for steps (i) and 
(iv) are shown in the left and right columns of Fig.0 
The contributions (ii) are found by calculating three- 
loop Green functions with operator insertions. Sample 
diagrams with a double insertion of dimension-six oper- 
ators are shown in the center column of Fig.|2| 

The Z-penguin contribution can be trivially obtained 
from that in K + — > n + vv, which has been recently com- 
puted at NNLO [13, E3- The electroweak box contri- 
bution on the other hand is slightly different for Kl — > 
and K + — * ■n Jr vv since the lepton line in the cor- 
responding Feynman diagrams is reversed and thus re- 
quires a new calculation. A comprehensive discussion of 
the technical details of the matching and the renormal- 
ization of the effective theory can be found in • 



TABLE II: The coefficients Kijkl arising in the approximate for- 
mula for P c at NNLO. 



^1000 


= -0.5373 


^0100 


= -6.0472 


^00 10 


= -0.0956 


KfJOOl 


= 0.0114 


KllOO 


= 3.9957 


K1010 


= 0.3604 


ftono 


= 0.0516 


«0101 


= -0.0658 


K2000 


= -0.1767 


^0200 


= 16.4465 


^0020 


= -0.1294 


^0030 


= 0.0725 



In the following we present only the final result for the 
0(a 2 ) matching correction C^ 2 \ the 0(a 3 s ) anomalous 
dimension yjt , and the 0(a 2 ) matrix element r^ 2 \ 
Employing the operator basis of 0, we obtain for the 
standard choices of Casimir operators Ca — 3, Cp = 4/3, 
and / active quark flavors 



-3(2) 



416 

~3~ 



16tt 2 



272 

3 M 2 



16 In 



2 

Ml 



Rm 27032 , x 1040 

7f ) = ^7-- 1088 C(3)-^/, 



9 



9 



(5) 



^ (2) = _112_80 ln ^_ i61n2 ^ 



Here C( x ) is the Riemann zeta function with the value 
£(3) w 1.20206 and m c = m c (/i c ) denotes the charm- 
quark MS mass. Our results for the NLO Wilson coeffi- 
cient, the anomalous dimension and the matrix element 
agree with the findings of |ll| where an error made in the 
original calculation [5j has been corrected. 

The analytic expression for P c including the complete 
NNLO corrections is too complicated and too long to 
be presented here. Instead setting A = 0.225, m t (m t ) = 
162.8 GeV and fi w = 80.0 GeV we derive an approximate 
formula for P c that summarizes the dominant parametric 
and theoretical uncertainties due to m c (m c ), a s (M z ), ji c , 
and /ib- It reads 



P r = 0.1198 



/ \ \ 2.3595 / ;„ n 6.6055 

m c (m c ) \ ( a s (M z ) 



1.30 GeV 



0.1187 



(6) 



1 + K ijhnLl na L J as L^ li: L l ^ b 



where 



hi 



In 



m c (m c ) 
1.30 GeV 

1.5 GeV 



1 0.1187 

L,,. = In ( — — z — 
Mb V 5.0 GeV 



(7) 



and the sum includes the expansion coefficients Kijki 
given in TablelTTl The above formula approximates the 
exact NNLO result with an accuracy of better than 
±1.0% in the ranges 1.15 GeV < m c (m c ) < 1.45 GeV, 
0.1150 < a s (M z ) < 0.1230, 1.0 GeV < \i c < 3.0 GeV 
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and 2.5 GeV < fib < 10.0 GeV. The uncertainties due 
to m t (m t ), fj, w and the different methods of computing 
a s (fi c ) from a s (M z ), which are not quantified above, are 
all below ±0.2%. Their actual size at NNLO will be dis- 
cussed below. 

Using the input parameters listed in TableQ] we find 
at the NNLO level 



P c = 0.115 ± 0.008 th cor ± 0.008 mc ± 0.001 Qs 



(0.115 ±0.018) 



0.225 
A 



(8) 



where now the residual scale ambiguities and the uncer- 
tainty due to m c (m c ) are of the same size. Comparing 
these numbers with Eq. © we observe that our NNLO 
calculation reduces the theoretical uncertainty by a fac- 
tor of more than 3. 

As can be nicely seen in the lower plot of Fig.^ P c 
depends very weakly on fi c at NNLO, varying by only 
±0.007. Furthermore, the three different treatments of 
a s affect the NNLO result in a negligible way. The three- 
loop values of a s (fi c ) used in the numerical analysis have 
been obtained with the program RUNDEC. The theoretical 
error quoted in Eq. JHJ includes also the dependence on 
fib and fi w of combined ±0.001. The presented scale 
uncertainties correspond to the ranges given earlier. 

Using Eqs. {Q, J3J), and JSJ the result in Eq. (@J is mod- 
ified to the NNLO value 

B(K L -> fi + fi-) SD = (0.79 ± 0.04 Pc ± 0.08 othcr ) x 10~ 9 
= (0.79 ±0.12) x 10" 9 . (9) 

Obviously, at present the errors from M t , m c (m c ) and 
the CKM parameters veil the benefit of the NNLO cal- 
culation of P c presented in this Letter. 

Provided both P c and B(Kl — > fi + p~)sr> are known 
with sufficient precision useful bounds on the Wolfenstein 
parameter p can be obtained In particular for the 
measured branching ratio B(Kl — > fi + P~)sd close to its 
SM predictions, one finds that given uncertainties cr(P c ) 
and <j(B(Kl — » fi + /j~)sd) translate into 



^1 = ± o. 89 ^) ±2.59^^ 



A^V )sd) 



B{K L -> fi+fi-) 



SD 



(10) 



As seen in Eq. (|10fl the accuracy of the determination 
of p depends sensitively on the error in P c . The reduction 
of the theoretical error in P c from ±22% down to ±7% 
translates into the following uncertainties 



cr(p) (±20%, NLO, 
~P~ ~ \±6% : 



NNLO. 



(11) 



implying a significant improvement of the NNLO over the 
NLO result. In obtaining these numbers we have included 
only the theoretical errors quoted in Eqs. © and |J8}. 



Using the conservative upper bound 

B{K L -» p + p-) SD < 2.5 x 10" 



(12) 



on the short-distance part of the Kl — > fi + fi branching 
ratio derived in Q , we find the following allowed range 



-0.74 < p < 3.13, 



(13) 



for the Wolfenstein parameter p employing a customized 
version of the CKMFITTER code. 

To conclude, we have evaluated the complete NNLO 
correction of the charm-quark contribution to B(Kl — * 
p + fi~)sr>- The inclusion of these contributions leads to 
a drastic reduction of the theoretical uncertainty in the 
relevant parameter P c . This strengthens the power of the 
rare decay Kl m determining the Wolfenstein 

parameter p from its short-distance branching ratio. 
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